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Theoretical calculations of neutralino-nucleon interaction rates with various nuclei are of great 
interest to direct dark matter searches such as CDMS, EDELWEISS, ZEPLIN, and other exper- 
iments since they are used to establish upper bounds on the WIMP-proton cross section. These 
interaction rates and cross sections are generally computed with standard, one or two parameter 
model-dependent nuclear form factors, which may not exactly mirror the actual form factor for 
the particular nucleus in question. As is well known, elastic electron scattering can allow for very 
precise determinations of nuclear form factors and hence nuclear charge densities for spherical or 
near-spherical nuclei. We use charge densities derived from elastic electron scattering data to cal- 
culate model independent, analytic form factors for various target nuclei important in dark matter 
searches, such as Si, Ge, S, Ca and others. We have found that for nuclear recoils in the range of 
1-100 keV significant differences in cross sections and rates exist when the model independent form 
factors are used: at 30 keV nuclear recoil the form factors squared differ by a factor of 1.06 for 28 Si, 
1.11 for 40 Ca, 1.27 for 70 Ge, and 1.92 for 129 Xe. We show the effect of different form factors on the 
upper limit on the WIMP-proton cross section obtained with a hypothetical 70 Ge detector during 
a 100 kg-day effective exposure. Helm form factors with various parameter choices differ at most 
by 10-20% from the best (Fourier Bessel) form factor, and can approach it to better than 1% if the 
parameters are chosen to mimic the actual nuclear density. 



Neutralinos, or more generically weakly-interacting massive particles (WIMPs), may be detected in the laboratory 
"directly" as they elastically collide with nuclei in a target, depositing enough energy to give typical nuclear recoils 
of keV. See [I|-[l3| for example for a review of calculating detection rates of neutralino dark matter. 

Direct detection experiments set limits on the WIMP-proton or neutralino-proton cross section in the following 
manner. First, the spin-independent neutralino-nucleus elastic cross section with a pointlikc nucleus of Z protons and 
N neutrons can be written as 



Hi — mM/(m + M) is the WIMP-nucleus reduced mass, with m the neutralino mass and M the nucleus mass. In 
most theoretical models it is assumed that G§ = G", so that the cross section scales with the square of the nucleus 
atomic number A. Thus the spin independent WIMP-proton cross section can be written as 



For a given detector, the expected number of recoil events with recoil energy in the range (E±,E2) is the sum over 
the nuclear species in the detector given by 
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where £i{E) is the effective exposure of each nuclear species in the detector (introduced in [I4I), and dRi/dE is the 
expected recoil rate per unit mass of species i per unit nucleus recoil energy and per unit time. The effective exposure 
is given by 

£ i =M i T i e i (E), (4) 

where 7$ is the active time of the detector during which a mass Aii of nuclear species i is exposed to the signal, and 
€i(E) is the counting efficiency for nuclear recoils of energy E. The differential rate dRi/dE is given by 

dR, _ paf 1 \F(q)\ 2 f f(v,t). 
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Here E is the energy of the recoiling nucleus, p is the local halo WIMP density, f(v,t) is the WIMP velocity 
distribution function in the frame of the detector, af 1 is the spin-independent WIMP-nucleus elastic cross section off 
a pointlike nucleus, and l-F^g)! 2 , with q = \J1ME (the recoil momentum), is a nuclear form factor. The upper bound 
on the WIMP-proton cross section is computed by comparing the expected number of events given by ([210 to the 
observationally set limit on the number of events detected in the relevant energy range(s). 

The determination of an upper limit on the WIMP-proton cross section depends of course on the shape of the signal 
in dR/dE, which in turn depends on two factors: 1) the integral over the WIMP velocity distribution, which contains 
the astrophysical uncertainties and which we will not discuss here, and 2) the form factor ^(g)! 2 , which contains 
the nuclear physics uncertainties and which we will examine. Notice that the upper limit set by direct detection 
experiments on <r p is inversely proportional to ^(q)! 2 . So a change of |F(g)| 2 by a factor of two implies an upper 
limit on a p stronger by a factor of two. It is thus crucial to understand and have a correct value for the nuclear form 
factor for the relevant range of momentum transfers. 

The nuclear form factor, F(q), is taken to be the Fourier transform of a spherically symmetric ground state mass 
distribution normalized so that F(0) = 1: 

F(V) = if f Pmass(r) e ~ tqr rf 3 r = i- f Pmass (r) ^r 2 dr. (6) 

M J M J qr 

Since the mass distribution in the nucleus is difficult to probe, it is generally assumed that mass and charge densities 
are proportional 

M 

PmassM = — PchargetV), (?) 

so that charge densities, determined through elastic electron scattering, can be utilized instead. Because of the 
normalization at q — 0, the proportionality assumption amounts to 

•Pmass('?) = ^chargc('Z)- (8) 

It is of course convenient to have an analytic expression for the form factor. Up to now, this expression has been 
provided by the Helm form factor [l8j given by 



where 
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is a spherical Bessel function of the first kind, and where R\ is an effective nuclear radius and s is the nuclear skin 
thickness, parameters that need to be fit separately for each nucleus. The Helm form factor was introduced as a 
modification of the form factor for a uniform sphere multiplied by a gaussian to account for the soft edge of the 
nucleus (see the next section for a more complete description) . The parameters R\ and s were in the past chosen to 
match numerical integration of Two-Parameter Fermi (Woods-Saxon) or other parametric models of nuclear density. 

In the past, Lewin and Smith [Io| demonstrated a method for fitting parameters in the Helm form factor to muon 
spectroscopy data in the Fricke et al. compilation [23| . They performed a two-parameter least squares fit to the muonic 
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spectroscopy data, finding the values of R\ and s for the Helm form factor which best reproduce the numerical fourier 
transform of a Two-Parameter Fermi distribution. Explicitly they set 



i?x = ^c 2 + 7 --k 2 cl 2 -bs 2 (11) 

and take s ~ 0.9 fin, a ~ 0.52 fm (as did Fricke et al. in their table IIIA), and 

c ~ 1.23^ 1/3 - 0.60 fin (12) 

(which is a least-square fit to the same table in Fricke et al.). This procedure, however, should be approached with 
caution due to the fact that the results depend on the nuclear density model (in this case Two-Parameter Fermi) 
which was used in the original fit to the data. Also, in the Fricke et al. compilation [23| in their table IIIA the value 
of the skin thickness was fixed, leading to a one parameter fit. Hence the fit to the form factor generated from the 
muon spectroscopy data is in essence a fit to a fit. The chief advantage here is that a more accurate Helm form factor 
is generated which is analytic and eliminates the need for numerical integration. 

It is true that the necessity of introducing a nuclear form factor correction, particularly in direct dark matter 
detection, has been widely recognized and written about in the literature (see for example [9j and [1Q(). However, 
what is not widely known in the dark matter community is that there exist model independent form factors derived 
from clastic electron scattering data that are both analytic (reproducing the chief advantage of using the Helm form 
factor) and more accurate than standard Helm form factors. These model independent form factors arc derived 
directly from elastic scattering data, and more importantly the relevant paramaterizations exist for the large range 
of nuclei relevant to dark matter searches (whereas fits to other paramaterizations such as the harmonic oscillator 
or modified harmonic oscillator model exist for only a few nuclei). We would like to suggest the use of these model 
independent, analytic expressions and point out that for many nuclei relevant to dark matter searches there exist 
significant differences between these form factors and the standard Helm parameterization. A new result which 
we will present is that these differences in form factors can lead to 10-20% shifts in the upper limits set on the 
WIMP-proton cross section as published by various experimental groups. 



II. NUCLEAR CHARGE DENSITIES AND FORM FACTORS 



Since the interaction between large nuclei and heavy neutralinos cannot be treated as scattering from point-like 
particles, it becomes necessary to introduce models for nuclear charge density into dark matter detection rate calcu- 
lations. In our review of the literature we found considerable difference in notation regarding various nuclear charge 
density/form factor models; we therefore begin with a brief review of the pertinent parameterizations. 

The simplest example of a model for the charge density of a nucleus is the uniform model in which the charge 
density is constant up to some cut-off radius R, i.e. 

Pair) = { ; < I (13) 

where the charge density is normalized such that the total charge contained in the nucleus is Ze. The form factor for 
this uniform model is simple to calculate and is given by 

Mi) = ^i(i R )- ( 14 ) 

Obviously such a charge density is nonphysical; a nucleus cannot exhibit such an infinitely sharp cutoff in its charge 
distribution. 

The Helm charge density solves the problem of the infinitely sharp cutoff in the uniform model by convoluting the 
uniform charge density with a Gaussian "surface smearing" density. The Helm charge distribution is given as 

Pn(r) = J ' Pu (r')p G (r-r ')&?', (15) 

where pg(^) is taken to be 



(16) 
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Here g is a parameter which is related to the radius of the Gaussian smearing surface. One advantage of the Helm 
charge density is that it has an extremely simple analytic form factor; in fact (by the convolution theorem) , the form 
factor is simply a product of the form factors of pu and pc [HI : 

F(q) = Fu(q)F G (q) 

= -^ji(qR)e-<?« 2 / 2 . (17) 

Returning to dark matter, in most dark matter direct detection calculations/simulations, the nuclear charge density 
is assumed instead to be of the Two-Parameter Fermi (Woods-Saxon) distribution form given by 

p(r) = -, — £ (18) 

' v 1 g(l — c)/o -L 1 V 1 

where c is the half-density radius, p c is the density at r = c, and the parameter a is related to the surface thickness 
t by t — (41n3)a; the charge density is again normalized by requiring the total charge contained in the nucleus to 
be Ze. The Two-Parameter Fermi (Woods-Saxon) distribution has been favored due to its relative simplicity (only 
two free parameters) as well as the fact that parameterizations for many nuclei have been determined from nuclear 
scattering experiments. However, the Fourier transform of the Two-Parameter Fermi or Woods-Saxon distribution 
must be calculated numerically since no closed form analytical Fourier transform exists. For computational simplicity, 
the spin independent form factor is usually taken to be of the Helm form as described above. As described in Q and 
(Tl | the Helm form factor is "virtually indistinguishable" from the numerical Woods-Saxon/Two-Parameter Fermi 
form factor. For example, DarkSUSY 4.1 [l4| takes the spin independent form factor to be of the Helm form ^ with 

Ri = \/R 2 - 5s 2 , (19) 



R = [0.91 (Af/GeV) 1/3 + 0.3] fm, (20) 

and 

s = 1 fm. (21) 
Here the expression for R is that of the nuclear radius in a common parametrization (l5l . flt| , while s is taken from 

The value of the momentum transfer at which the form factor is evaluated depends on the recoil energy and the 
mass of the nucleus in question. The momentum transfer can be written as 



q = V2ME, (22) 

where M is the mass of the target nucleus and E is the nuclear recoil energy. The momentum transfer q, commonly 
given in GeV, can easily be converted to fm -1 through 



! _ y/2M[GeV]£[keV] x 10-6 
q[ J ~ 0.197 GcVfm ' [26) 



where the nuclear mass is given in GeV and the nuclear recoil energy in keV. 



III. ELASTIC ELECTRON SCATTERING AND MODEL INDEPENDENT FORM FACTORS 

Elastic electron scattering as a tool for probing nuclear structure was pioneered in 1953 at the Stanford Linear 
Accelerator by Hofstadter and collaborators [19j . In the next twenty to thirty years a tremendous amount of data 
was gathered which allowed precise determination of nuclear charge distribution parameters culminating in a 1986 
compilation of nuclear charge density parameters for over one hundred nuclei in Atomic Data and Nuclear Data 
Tables pol ]. Many nuclei have charge densities expressed in multiple parameterizations: two parameter Fermi, three 
parameter Fermi, three-parameter gaussian, etc. 

Although the Two-Parameter Fermi or Woods-Saxon charge distribution fits many nuclei well, actual nuclear charge 
densities may be more complex than this simple parameterization allows for. For example, Figure 1 shows the charge 
densities for 40 Ca and 208 Pb plotted versus nuclear radius. 40 Ca shows an increase in density towards the nuclear 
center, while 208 Pb shows a non-constant interior density manifesting in a depressed charge density at about 3 fm (see 
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FIG. 1: Charge density for Ca and Pb as obtained from elastic electron scattering data. 



for example [25| and [26[). These features are not well-produced by a generic Two-Parameter Fermi or Woods-Saxon 
charge distribution. More complex charge densities are therefore necessary to reproduce these features. 

Additionally, most charge densities determined from electron scattering experiments have been analyzed using 
model distributions such as the Two-Parameter Fermi or Woods-Saxon distributions. The drawback here is that 
charge density parameters measured are fundamentally model dependent. This point applies as well to harmonic 
oscillator or modified harmonic oscillator charge density parameterizations, despite the fact that these have been used 
to characterize a limited number of nuclei involved in dark matter searches. As noted by Sick [26|, it is at times 
unclear whether the charge density parameters measured are the fundamental characterization of the nuclear charge 
density or if they are simply the result of fitting data to a model which may be too restrictive in form. To combat this 
fundamental uncertainty, model independent methods for analyzing electron scattering data and extracting nuclear 
charge densities were developed by several groups. Of course, these methods were not completely model independent; 
some model dependence is necessary since electron scattering data is only available over a finite momentum transfer 
range. Two primary methods emerged for extracting nuclear charge density parameters in a model independent 
fashion: an approach in which the charge density is written as a sum of gaussians and an alternate approach in which 
the charge density is written as a Fourier superposition of Bessel functions. 
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In the Sum of Gaussians expansion (SOG), first introduced by Sick [26|], the charge density of a nucleus is modelled 
as a series of Gaussians. The width of the Gaussians, given by the parameter 7, is calculated by equating it to the 
smallest width of the peaks in the nuclear wave functions as found in Hartree-Fock calculations; negative amplitudes 
for the Gaussians are not allowed as to avoid creating structures with widths smaller than 7 due to interference. As 
noted in [20J this method has several advantages: values of p(r) at different radii are decoupled due to the rapid 
fall-off of the Gaussian tail, and as long as a sufficient number of Gaussians are used to ensure a good fit, the results 
of the analysis is independent of the number of Gaussians. 

In the SOG expansion the charge density is written as 



p(r) = f> {e-K-* 



)hY +e -[(r+Ri)hY 



(24) 



where the coefficients Ai are given by 



A, 



ZeQt 



2 7 r 3 / 2 7 3 (l + 2 J Rf/7 2 )' 



(25) 



The Qi represent the fractional charge carried in the i th gaussian and lead to the following definition for the normal- 
ization of the charge density: 



A' 



(26) 



An analytical form factor can be determined for this density parameterization, which eliminates the necessity of 
performing numerical integration to find the form factor as in the Two-Parameter Fermi density parameterization. 
Assuming spherical symmetry, the form factor in the SOG expansion is given by 



F(q) 



N 

V/*r_ 
»=i 



Qi 



2i?f/ 7 2 



cos(qRi 



2Rj sin(qRi) 



r 



qRi 



(27) 



The other model independent parameterization of nuclear charge densities that was developed is a Fourier-Bessel 
expansion. In the Fourier-Bessel expansion (FB), first introduced by Dreher et al. (27[, the charge density is modelled 
as a sum of Bessel functions up to some cut-off radius i?, and is assumed to be zero thereafter. The form factor is 
assumed to fall off for large q as q~ A and e~ aq ; these result from the distribution of nucleons in the nucleus and from 
the finite extension of the nucleons respectively. Although the results depend slightly on the value of R, the cut-off 
radius, an advantage of this approach is that this collection of assumptions gives an upper limit on the contribution 
of higher Fourier components to the series expansion [2l| . 

Specifically, 



p(r) 



Y^=i Qujo (vwr/R) for r < R, 
for r > R, 



(28) 



where jo(x) = smx/x is the zero-th order spherical Bessel function. Here the charge density is normalized by requiring 



p(r)d 3 r = Ze. 



(29) 



This charge density also has the advantage that an analytical expression may be found for the form factor. Assuming 
spherical symmetry, Fourier transforming the charge density yields 



N 

E 



F(q) 



sin(gfl) gj ^ - 1 2R2 
qR " (-!)»«,, 

i/=l 



(30) 



which is normalized to F(0) = 1. 

In a final point regarding model independent charge densities, it is difficult to estimate the theoretical errors 
involved in fits to elastic electron scattering data. As noted in [20| . in the case of the FB approach the errors in 
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the individual coefficients are not presented in data compilations since errors are strongly correlated and can only be 
extracted from the full correlation matrix; as this matrix is never published in papers, it is impossible to give errors on 
individual coefficients of the FB expansion for the charge density. Recently, however, Anni, Co', and Pellegrino in [28[ 
analyzed model independent extraction of charge density parameters from elastic electron scattering data to determine 
uncertainties and the minimum number of expansion coefficients needed to give an accurate representation of the data. 
Anni et al. determined that in the model independent approach a truncation error is unavoidable; however, they 
developed an approach to determine an optimal number of coefficients. They point out that distributions extracted 
from these models must be used with care; however, these models are still much more robust than model dependent 
fits such as the Woods-Saxon or Two-parameter Fermi model. 



IV. NUCLEAR CHARGE DENSITY PARAMETERS FROM MUONIC ATOM SPECTROSCOPY 

One can also extract nuclear charge parameter factors from nuclei using muonic atom spectroscopy, however pa- 
rameters extracted in this manner are to some degree model dependent as the final parameter fitting depends on the 
choice of an analytical charge density. The nuclear charge distributions are extracted by considering finite size effects 
to the energy shift in first-order perturbation theory (see [22| and (23|, for example, for a review). The analysis takes 
advantage of Barret moments [24j which can be extracted in a model independent fashion from the transition energies. 
Charge density parameters are determined by finding the eigenvalues of the Dirac equation with the analytic charge 
density fitted to the experimental transition densities (23[. In the analysis of muonic atom spectroscopy data the 
Two-Parameter Fermi charge distribution is used: 

P( r ) = -T—T7 7- (31) 

The skin thickness t = 4aln(3) is fixed to 2.30 fm, and the half-density radius parameter c is fitted to reproduce 
the experimental transition energies. This method can also be applied to deformed nuclei by writing the half-density 
radius c as 

c = Ro[l + lhY m (e,4>)], (32) 

where fa is the quadrupole deformation parameter, and Rq is the monopole radius. 

As mentioned previously, the disadvantage to the Two-Parameter Fermi charge density lies in the fact that it has 
no analytic Fourier transform, and thus form factors must be computed through numerical integration. However, as 
noted in , data from muonic spectroscopy parameterizing the nuclear charge density in terms of a Two-Parameter 
Fermi distribution exists for several nuclei relevant to direct dark matter detection experiments such as Na, Xe, 
and I. Two-Parameter Fermi charge distribution parameters also exist for 184 > 186 W. However, these parameters were 
obtained from elastic electron scattering experiments. In order to determine the impact of using the Helm versus the 
form factor associated with the Two-Parameter Fermi charge distribution (and to compare against utilization of the 
Sum of Gaussian or Fourier-Bessel parameterizations) we numerically fourier transform the Two-Parameter Fermi 
distribution obtaining the associated form factor, generically referring to it as a Woods-Saxon form factor. Numerical 
integration was performed with a 48-point Legendre-Gauss quadrature. Although integrating the oscillating sin(qr) 
in [5] is notoriously subtle, the 48-point Legendre-Gauss qaudrature used gave robust results and reproduced form 
factors previously published in the literature (for example, see 10]). The half-density parameter c, the parameter a 
(which is related to the nuclear skin thickness t) , and the rms value of the charge density for the nuclei of interest are 
included the the appendix as Table 7. 



V. DISCUSSION AND RESULTS 



To determine the effect of using model independent form factors on dark matter direct detection rates, DarkSUSY 
was modified to compute Helm form factors using the standard parameterizations from [H[-[l3] (referred to here 
onward as Helm/DarkSUSY-4.1) and the Helm form factors of Lewin and Smith (referred to here onward as 
Helm/Lewin-Smith), as well as model independent form factors derived from the SOG or FB approach. The form 
factors were evaluated at the appropriate momentum transfers for nuclear recoils between 10 and 100 keV. The 1986 
compilation of electron scattering data [20j was used to calculate the relevant charge densities and form factors. Nuclei 
were chosen based on their importance to current and future/planned dark matter direct detection searches as well 
as the availability of electron scattering data in the appropriate momentum transfer range. We have restricted our 
attention to spin independent scattering of neutralinos from nuclei. This is not too restrictive since many of the nuclei 
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FIG. 2: Helm, FB, and SOG form factors for 28 Si versus nuclear radius as obtained from elastic electron scattering data. Lines 
for momentum transfers corresponding to 10 and 100 keV nuclear recoil are included. 



relevant in direct dark matter searches are even-even nuclei for whom the spin dependent cross section vanishes in 
the J = ground state (see for example Q). Since the relevant nuclei tend to be spherical or roughly spherical in 
shape, their charge densities and form factors are able to be analyzed using the SOG or FB approach. Iodine and 
Sodium (important to the DAMA experiment; see for example [111), Xenon (important to the ZEPLIN and XENON 
experiments; see for example [30l | and (3lj|). and Tungsten (important in the CRESST experiment; see for example 
[3a ]) are analyzed using the Two-Parameter Fermi charge density as in [l(| for the sake of completeness. 

A range of nuclei from 12 C to 208 Pb were analyzed to give a general trend for increasing nuclear mass. Of particular 
importance to dark matter experiments are the following nuclei: 28 Si (used in the CDMS experiments [HI), 32 S (used 
in the DRIFT experiment 34]), 40 Ca (used in the CRESST experiments [H), and 7 °- 74 Ge (used in the CDMS, 
EDELWEISS, GENIUS, and CryoArray experiments [32|-[36]). For a convenient summary chart of current and 
future-planned dark matter direct detection experiments and the associated target materials see Tables 1-3 in [371 ]. 

In Figure[2]we plot |-F(g)| 2 for 28 Si, an important target medium in the CDMS experiments, for small momentum 
transfers less than 1 fm _1 . Lines indicate the momentum transfers which corresponds to 10 keV and 100 keV 
nuclear recoils. As can be seen from the plot, all of the form factors in their various parameterizations are not 
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FIG. 3: Charge density for Ge in the Woods-Saxon and FB parameterizations. 



significantly different at low (~10 keV) nuclear recoils (they differ by at most 2%), and the Helm form factor (in 
either parameterization) may be used with confidence. However, for larger momentum transfers the FB, SOG, and 
Woods-Saxon form factors begin to diverge from the Helm form factors; at a momentum transfer corresponding to 
100 keV there is a 20.9% difference between the Helm/DarkSUSY-4.1 and the FB/SOG form factors. The Helm 
form factor of Lewin and Smith differs from the FB or SOG form factors negligibly at small momentum transfers 
and up to 5% at 100 keV nuclear recoil. Hence for small nuclei at low momentum transfers all parameterizations 
of the form factor fit well. One should also note, however, that the FB, SOG, and WS form factors, despite the 
different parameterizations, are virtually indistinguishable, giving further confidence in the use of model independent 
form factors. The Woods-Saxon form factor is a numerical fourier transform of a physical charge density model (the 
Two-Parameter Fermi), and should be in some sense thought of as the closest approximation to the actual form factor. 
The two model independent form factors (FB and SOG) trace the nuclear density more accurately, and easier to use 
as they are analytic: they should be preferred whenever available. 

In Figure [3] we plot the charge density for 70 Ge using both model independent fits to elastic electron scattering 
data (the FB fit) and also using the standard Woods-Saxon (Two-Parameter Fermi) charge density. As in Figure [IJ 
the charge density for germanium shows a non-constant interior density which cannot be characterized entirely by the 
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Woods-Saxon distribution. Figured] shows the Helm/DarkSUSY-4.1, Helm/Lewin-Smith, Woods-Saxon, and FB form 
factors for 70 Ge plotted over a range of to 4 fm _1 . As can be seen from the plot, the Helm/DarkSUSY-4.1 and the 
other form factors differ significantly over the range of momentum transfers. Most significantly, the first diffraction 
minimum in the FB and WS form factors occurs close to 1 fm _1 whereas the first minimum in the Helm/DarkSUSY- 
4.1 form factor occurs at significantly larger momentum transfer. The difference in form factors can be traced to the 
use of the nuclear radius formula (f2Uf , from [15| , instead of the more accurate one in Eq. (fTTj) , from [l(| ■ 

Figure [5] shows the same plot of form factors for 70 Ge but this time highlighting the relevant range of momentum 
transfers important for direct detection experiments. The momentum transfers corresponding to 10 and 100 keV 
nuclear recoils are highlighted by vertical lines. Although the divergence between the two form factors at 10 keV nuclear 
recoil is small (8.12% exactly between Helm/DarkSUSY-4.1 and FB form factors and 2.0% for the Helm/Lewin-Smith 
and the FB form factors), the form factors are already beginning to diverge. For larger nuclear recoils the form factors 
continue to diverge leading to a 27.4% and 6.3% difference at 30 keV, a 65.8% and 12.1% difference at 60 keV, and a 
148.0% and 19.6% difference at 100 keV for the Helm/DarkSUSY-4.1 vs. the FB and the Helm/Lewin-Smith vs. the 
FB form factor respectively. Again the WS and FB form factors are indistinguishable over the relevant momentum 
transfer range. For illustrative purposes, we plot the ratio of the form factor squared in the Helm scheme to that of 
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the form factor squared in the FB expansion for 70 Ge in Figure [6l 

In Figure [7]we plot |_F((7)| 2 for 208 Pb for momentum transfers in the range of zero to about 2 fm _1 . Although no 
direct dark matter searches currently employ lead as a detection medium (though it is used for shielding), lead is an 
excellent example of the general trend for heavier nuclei. As the atomic number and mass of the nucleus increases, 
the Helm/DarkSUSY-4.1, FB, SOG, and WS form factors begin to diverge at smaller momentum transfers. As can 
be seen in the case of 208 Pb, the first diffraction minimum for the FB form factor occurs between 0.6 and 0.7 fin -1 
whereas the first diffraction minimum for the Hclm/DarkSUSY-4.1 form factor occurs for the much larger momentum 
transfer of about 0.9 fm _1 . The Helm form factor of Lewin and Smith does a much better job matching the WS, FB, 
and SOG form factors, however the percent difference can become quite large (as high as 41% as 60 keV recoil energy). 
This discrepancy leads to a large difference in the spin independent cross sections of neutralinos with 208 Pb for nuclear 
recoils as small as 10 keV. Although these diffraction minima in electron or muon scattering can be partially filled by 
multiple photon exchanges in the nucleus, in the case of neutralinos this is not expected to occur [lOj. Lead clearly 
illustrates the hazards of using large-A nuclei as a detection medium; potential signals may be cut-off due to very 
small form factors in the relevant range of momentum transfers. 

In summary, Table 1 gives the ratios the Helm form factor squared to either FB, SOG, or WS form factors squared 
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FIG. 6: Ratio of the form factors squared in the helm and Fourier Bessel schemes for 70 Ge. 



for various nuclei important in dark matter searches at 10, 30, 60, and 100 keV nuclear recoils. Note that 27 A1 
is unique in that its Helm form factor as computed in DarkSUSY 4.1 is actually smaller than the corresponding 
model independent form factor; expected detection rates for an Al target should therefore be higher than previously 
expected. Table 2 gives the ratio of the Hclm/Lewin-Smith form factor squared to either the FB, SOG, or WS form 
factors squared for various nuclei at a series of nuclear recoil energies. 

In Table [H the starred entries show an extreme difference between the Helm and other form factors in which the 
ratio is greater than one hundred; this occurs when the form factor is evaluated on or about a diffraction minimum 
for the WS, FB, or SOG form factors with the corresponding Helm/DarkSUSY-4.1 form factor diffraction minimum 
occurring at higher momentum transfer. Note that no such entries occur in Table 2 as the Helm/Lewin-Smith form 
factor has diffraction minima whose placement roughly matches the locations in the FB, SOG, or WS form factors. 
As can be seen from Tables Q] and HH although either Helm form factor may be used with confidence for light nuclei 
(A < 30) at low nuclear recoils, model independent form factors become increasingly more accurate for large nuclei 
even at relatively modest momentum transfers. Failure to account for this correction to the form factor can lead to 
10-20% errors in cross sections and detection rates. 

We want to illustrate the importance of using accurate form factors when setting limits on WIMP-proton cross 
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FIG. 7: Helm and FB form factors for 208 Pb with 10 and 100 keV recoil lines for momentum transfers less than 3.0 fm 

sections. For this purpose, we imagine a hypothetical detector made of 70 Ge. We assume a simple isothermal sphere 
model for the WIMPs in the Milky Way halo, so that the WIMP speeds with respect to the halo obey a Maxwellian 
distribution with a velocity dispersion truncated at the escape velocity. We use the same halo velocity distribution 
function as [38[ - see particularly (30)-(33) - setting ah = 220 km/s, ve = 220 km/s, v csc = 650 km/s, and p = 0.3 
GeV/c 2 /cm 3 . We assume a detector energy threshold of 10 keV, an efficiency of 50%, an exposure of 100 kg-days, 
and no events observed in the energy range of 10-50 keV. 

In Figure [5] we show the upper limit on the WIMP-proton cross section obtained using (J2H5]) with A^o-sokcV < 2.3 
and various models of the form factor. The form factor as given by DarkSUSY 4.1 is a Helm type form factor 
calculated with the parameters given in (|19ti21|) . The form factor computed using the procedure in Lewin and Smith 
[l(| is a Helm type form factor with parameters given by (|1HI12[) . The Fourier Bessel form factor is calculated as in 
Equation [3D] using the data given in Appendix I Table 4. Finally, the form factor labelled "Helm/2pF" is a Helm type 
form factor computed using R\ in (|11|) with s = 0.9 fm and with a = 0.5807 fm and c = 4.430 fm directly from the 
Two-Parameter Fermi distribution in table VIII in Fricke et al. (instead of the least-square fit to their table IIIA in 
Lewin and Smith). 

In this example the limits computed with DarkSUSY 4.1 differ from the Fourier Bessel limits by 10-20%, while 
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Nucleus 


F 2 (q) ii /F 2 ( q ) 
10 keV 


F 2 (q)u/F 2 (q) 
30 keV 


F 2 (q)u/F 2 (q) 
60 keV 


F 2 (q)n/F 2 (q) 
100 keV 


12 C 


1.01 


1.01 


1.03 


1.06 


ICQ 


1.01 


1.03 


1.06 


1.10 


23 Na 


1.02 


1.05 


1.10 


1.17 


27 Al 


0.99 


0.97 


0.94 


0.90 


28 Si 


1.02 


1.06 


1.12 


1.21 




1.02 


1.08 


1.16 


1.28 


40 Ar 


1.03 


1.10 


1.22 


1.41 


40 Ca 


1.03 


1.11 


1.23 


1.42 


70 Ge 


1.08 


1.27 


1.66 


2.48 


127j 


1.21 


1.89 


4.98 


13739.7* 


129 Xe 


1.22 


1.92 


5.36 


1164.1* 


134 Xe 


1.22 


1.96 


5.91 


100.8* 


184 W 


1.48 


4.68 


20.80 


0.02 


208p b 


1.54 


6.91 


1.64 


1.15 



TABLE I: Ratios of the Helm/DarkSUSY 4.1 [TiJ to FB or SOG or WS form factors squared for various nuclei important in 
current or future direct dark matter searches at 10, 30, 60 and 100 keV nuclear recoil. Lead, though not important in direct 
dark matter searches, is added as an example of a very large A nucleus. 



Nucleus 


F 2 (q) ll /F 2 (q) 
10 keV 


F 2 (q) K /F 2 (q) 
30 keV 


F 2 (q) H /F 2 (q) 
60 keV 


F 2 (q) H /F 2 (q) 
100 keV 


12 C 


1.00 


1.00 


1.01 


1.01 


16 Q 


1.00 


1.01 


1.02 


1.03 


23 Na 


1.01 


1.02 


1.03 


1.05 


27 Al 


1.00 


1.01 


1.02 


1.04 


28 Si 


1.00 


1.01 


1.03 


1.05 


32 S 


1.01 


1.02 


1.05 


1.08 


40 Ar 


1.01 


1.03 


1.05 


1.09 


40 Ca 


1.01 


1.03 


1.06 


1.1 


70 Ge 


1.02 


1.06 


1.12 


1.19 


127j 


1.03 


1.08 


1.13 


1.27 


129 Xe 


1.03 


1.08 


1.13 


2.98 


134 Xe 


1.02 


1.05 


1.03 


2.67 




1.07 


1.25 


0.99 


1.6 


208p b 


1.03 


1.01 


1.41 


1.23 



TABLE II: Ratios of the Helm/Lewin- Smith form factor [l(| squared to FB or SOG or WS form factors squared for various 
nuclei important in current or future direct dark matter searches at 10, 30, 60 and 100 keV nuclear recoil. Lead, though not 
important in direct dark matter searches, is added as an example of a very large A nucleus. 

those computed a la Lewin and Smith differ from the Fourier Bessel limits by 1-2%. (Notice that the logarithmic 
vertical scale covers a range much smaller than the usual one in this kind of plots.) Numerically, we can compare 
our hypothetical upper limit on the WIMP-proton cross section for a 100 GcV WIMP, say. Using the Fourier Bessel 
form factor we find the limit to be 2.38 x 10 -7 pb. Using the Helm form factor with a and c parameters directly from 
Frickc et al. we find it to be 2.39 x 10~ 7 pb, which is 0.5% larger. Using the Helm form factor a la Lewin and Smith 
we find the limit to be 2.34 x 10 -7 pb, which is 1.5% smaller. Using the Helm form factor computed with DarkSUSY 
4.1 we find the upper limit to be 1.93 x 10 -7 pb, which is 19% smaller. 

The improvement in form factors that incorporate nuclear data more and more accurately is quite evident in this 
example. The upper limits on the WIMP-proton cross section calculated using form factors generated with the Helm 
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FIG. 8: Hypothetical upper limit on the WIMP-proton cross section for a 100 kg-day exposure of 70 Ge (assuming an efficiency 
of 50% and no detected event) and various form factors. 



distribution utilizing muonic spectroscopy data compiled in Fricke et al. [231 ] differ only by less than 1% from the 
WIMP-proton cross section calculated using the Fourier Bessel form factor. The Helm form factor computed a la 
Lewin and Smith is also a very good approximation to the Fourier Bessel form factor, to within a few percent. The 
Helm form factor in DarkSUSY 4.1 is further off, by 10-20%, because of a poor choice of the parameter Ri, which is 
there set equal to the nuclear radius in standard parameterizations (Eq. [2Ql from (Til]). 



VI. CHARGE VS. MASS FORM FACTORS 



Neutralinos, of course, interact with the mass distribution of the nucleus, and for neutralino scattering we should 
be using mass rather than charge form factors. Throughout this paper we have assumed (as have previous authors, 
for example Q) that the mass and charge distributions within nuclei arc roughly proportional. However, there have 
been recent attempts to pin down the nuclear mass distribution using coherent photoproduction of n° mesons (40j . 
Photoproduction of pions is in fact an attractive method for studying nuclear mass distributions; unlike induced 
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hadron reactions it probes the entire nuclear volume. In a recent paper on pion photoproduction [4l| Krusche 
analyzed the mass form factors of 12 C, 40 Ca and ™ at Pb using the Helm model for the form factor. He found that 
rms mass radii were slightly smaller than rms charge radii, though it is still unclear whether these results are real or 
merely poorly understood systematic effects in the data or in the model dependent calculations. As measurements of 
coherent photoproduction of pions from nuclei improve, more stringent limits of differences between the nuclear mass 
and charge densities will be set, including the ability to analyze the data in a model independent fashion as with 
electron scattering data. 

VII. CONCLUSION 

In conclusion, significant differences may exist even at relatively low momentum transfers between generic Helm- 
type and more realistic model independent form factors extracted from elastic electron scattering data, particularly 
for large A nuclei. These new form factors should be utilized to avoid errors in published limits on neutralino-nucleon 
cross sections and detection rates in direct dark matter searches. The Fourier Bessel and Sum of Gaussian form 
factors have the advantage that they are model independent, analytic, and are derived from electron scattering data; 
we have shown that for calculating upper bounds on the WIMP-proton cross section they are at least as accurate 
as modified Helm form factors which have been fit for a specific nucleus to muonic spectroscopy data. The FB 
and SOG form factors have the advantage of existing for the wide array of nuclei being used in current and future 
planned dark matter detectors, and these parameterizations can be used without the need for extensive fitting to 
experimental data for different nuclei. We suggest that their use is both simpler and more convenient than modified 
Helm form factors. Of course, at present, one still must assume that the charge distribution in the nucleus mirrors 
the mass distribution; however, photo-pion production experiments are beginning to probe nuclear interiors to give 
model independent parameterizations of the nuclear mass density and should provide a wealth of new data within the 
next few years. 

The model independent analytic form factors in both the Sum of Gaussians (SOG) and Fourier-Bessel (FB) ap- 
proaches have been incorporated into the DarkSUSY [14j code along with a numerical integration routine to calculate 
Fourier transforms of Two-Parameter Fermi distributions; these improvements should appear in the next major public 
release of the program and will replace the need to rely on the simpler Helm form factor. Parameters in the SOG and 
FB expansions for the most commonly used target nuclei are included in the appendix. 
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IX. APPENDIX I: NUCLEAR CHARGE DENSITY PARAMETERS FOR NUCLEAR FORM FACTORS 



Nucleus 


12q 


16q 


28 Si 


30 Si 


rms [fm] 


2.464(12) 


2.737(8) 


3.085(17) 


3.173(25) 


ai 


0.15721e-l 


0.20238e-l 


0.33495e-l 


0.28397e-l 


ai 


0.38897e-l 


0.44793e-l 


0.59533e-l 


0.54163e-l 




0.37085e-l 


0.33533e-l 


0.20979e-l 


0.25167e-l 


ai 


0.14795e-l 


0.35030e-2 


-0.16900e-l 


-0.12858e-l 


as 


-0.44831e-2 


-0.12293e-l 


-0.14998e-l 


-0.17592e-l 


aa 


-0.10057e-l 


-0.10329e-l 


-0.93248e-3 


-0.46722e-2 


07 


-0.68696e-2 


-0.34036e-2 


0.33266e-2 


0.24804e-2 


as 


-0.28813e-2 


-0.41627e-3 


0.59244e-3 


0.14760e-2 


ag 


-0.77229e-3 


-0.94435e-3 


-0.40013e-3 


-0.30168c-3 


aio 


0.66908e-4 


-0.25771e-3 


0.12242e-3 


0.483464e-4 


an 


0.10636e-3 


0.23759e-3 


-0.12994e-4 


O.OOOOOeO 


012 


-0.36864e-4 


-0.10603e-3 


-0.92784e-5 


-0.51570e-5 


ai3 


-0.50135e-5 


0.41480e-4 


0.72595e-5 


0.30261e-5 


ai4 


0.94550e-5 




-0.42096e-5 




ais 


-0.47687e-5 








7? [fm] 


8.0 


8.0 


8.0 


8.5 



TABLE III: Fourier-Bessel coefficients for 12 C, le O, 28 Si and 30 Si as in [IJ, to be used in (j30j. 



Nucleus 


32g 


40 Ar 


40 Ca 


70 Ge 


rms [fm] 


3.248(4) 


3.423(14) 


3.450(10) 


4.043(2) 


01 


0.37251e-l 


0.30451e-l 


0.44846e-l 


0.38182e-l 


a 2 


0.60248e-l 


0.55337e-l 


0.61326e-l 


0.60306e-l 


a-3 


0.14748e-l 


0.20203e-l 


-0.16818e-2 


0.64346e-2 




-0.18352e-l 


-0.16765e-l 


-0.26217c-l 


-0.29427e-l 


as 


-0.10347e-l 


-0.13578e-l 


-0.29725e-2 


-0.95888e-2 


ae 


0.30461e-2 


-0.43204e-4 


0.85534e-2 


0.87849e-2 


a 7 


0.35277e-2 


0.91988e-3 


0.35322e-2 


0.49187e-2 


as 


-0.39834e-4 


-0.41205e-3 


-0.48258e-3 


-0.15189e-2 


ag 


-0.97177e-4 


0.11971e-3 


-0.39346e-3 


-0.17385e-2 


010 


0.92279e-4 


-0.19801e-4 


0.20338e-3 


-0.16794e-3 


an 


-0.51931e-4 


-0.43204e-5 


0.25461e-4 


-0.11746e-3 


012 


0.22958e-4 


0.61205e-5 


-0.17794e-4 


0.65768e-4 


«13 


-0.86609e-5 


-0.37803e-5 


0.67394e-5 


-0.30691e-4 


014 


0.28879e-5 


0.18001e-5 


-0.21033e-5 


0.13051e-5 


«15 


-0.86632e-6 


-0.77407e-6 




-0.52251e-5 


R [fm] 


8.0 


9.0 


8.0 


10.0 



TABLE IV: Fourier-Bessel coefficients for 32 S, 40 Ar, 40 Ca, and 70 Ge as in [HJ, to be used in (J30J). 
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Nucleus 


72 Ge 


74 Ge 


76 Ge 


208p b 


rms [fm] 


4.060(2) 


4.075(2) 


4.081(2) 


5.499(1) 


ai 


0.38083e-l 


0.37989e-l 


0.37951e-l 


0.62732e-l 


ai 


0.59342e-l 


0.58298e-l 


0.57876e-l 


0.38542e-l 


a-i 


0.47718e-2 


0.27406e-2 


0.15303e-2 


-0.55105e-l 


04 


-0.29953e-l 


-0.30666e-l 


-0.31822e-l 


-0.26990e-2 


a 5 


-0.88476e-2 


-0.81505e-2 


-0.76875e-2 


0.31016e-l 


06 


0.96205e-2 


0.10231e-l 


0.11237e-l 


-0.99486e-2 


a 7 


0.47901e-2 


0.49382e-2 


0.50780e-2 


-0.93012e-2 


08 


-0.16869e-2 


-0.16270e-2 


-0.17293e-2 


0.76653e-2 


ag 


-0.15406e-2 


-0.13937e-2 


-0.15523e-2 


0.20886e-2 


oio 


-0.97230e-4 


0.15376e-3 


0.72439e-4 


-0.17840e-2 


an 


-0.47640e-4 


0.14396e-3 


0.16560e-3 


0.74876e-4 


012 


-0.15669e-5 


-0.73075e-4 


-0.86631e-4 


0.32278e-3 


«13 


0.67076e-5 


0.31998e-4 


0.39159e-4 


-0.11353e-3 


014 


-0.44500e-5 


-0.12822e-4 


-0.16259e-4 




Ol5 


0.22158e-5 


0.48406e-5 


0.63681e-5 




i? [fm] 


10.0 


10.0 


10.0 


11.0 



TABLE V: Fourier-Bessel coefficients for 72 Ge, 74 Ge, 76 Ge, and 208 Pb as in [Ij, to be used in J30). 





12 C 


16 


28 Si 


rms [fm] 


2.469(6) 


2.711 


3.121 


i 


Ri 


Qt 




Qr 


Ri 


Qi 


1 


0.0 


0.016690 


0.4 


0.057056 


0.4 


0.033149 


2 


0.4 


0.050325 


1.1 


0.195701 


1.0 


0.106452 


3 


1.0 


0.138621 


1.9 


0.311188 


1.9 


0.206866 


4 


1.3 


0.180515 


2.2 


0.224321 


2.4 


0.286391 


5 


1.7 


0.219097 


2.7 


0.059946 


3.2 


0.250448 


6 


2.3 


0.278416 


3.3 


0.135714 


3.6 


0.056944 


7 


2.7 


0.058779 


4.1 


0.000024 


4.1 


0.016829 


8 


3.5 


0.057817 


4.6 


0.013961 


4.6 


0.039630 


9 


4.3 


0.007739 


5.3 


0.000007 


5.1 


0.000002 


10 


5.4 


0.02001 


5.6 


0.000002 


5.5 


0.000938 


11 


6.7 


0.00007 


5.9 


0.002096 


6.0 


0.000002 


12 






6.4 


0.000002 


6.9 


0.002366 


RP [fm] 


1.20 


1.30 


1.30 



TABLE VI: Sum of Gaussian coefficients for 12 C, 16 0, and 28 Si as in [2(J, to be used in ((27)) . RP, the rms radius of the 
Gaussians, is related to 7 by RP = 7^/3/2. 
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32 g 


40 Ca 


208p b 


rms [fm] 


3.258 


3.480(3) 


5.503(2) 


i 


Ri 


Qi 


Rr 


Qr 


R, 


Qi 


1 


0.4 


0.045356 


0.4 


0.042870 


0.1 


0.003845 


2 


1.1 


0.067478 


1.2 


0.056020 


0.7 


0.009724 


3 


1.7 


0.172560 


1.8 


0.167853 


1.6 


0.033093 


4 


2.5 


0.324870 


2.7 


0.317962 


2.1 


0.000120 


5 


3.2 


0.254889 


3.2 


0.155450 


2.7 


0.083107 


6 


4.0 


0.101799 


3.6 


0.161897 


3.5 


0.080869 


7 


4.6 


0.022166 


4.3 


0.053763 


4.2 


0.139957 


8 


5.0 


0.002081 


4.6 


0.032612 


5.1 


0.260892 


9 


5.5 


0.005616 


5.4 


0.004803 


6.0 


0.336013 


10 


6.3 


0.000020 


6.3 


0.004541 


6.6 


0.033637 


11 


7.3 


0.000020 


6.6 


0.000015 


7.6 


0.018729 


12 


7.7 


0.003219 


8.1 


0.002218 


8.7 


0.000020 


RP [fm] 


1.35 


1.45 


1.70 



TABLE VII: Sum of Gaussian coefficients for 32 S, 40 Ca, and 208 Pb as in [Hj, to be used in J27J}. RP, the rms radius of the 
Gaussians, is related to 7 by RP = 7^/3/2. 



Nucleus 


c [fm] 




a [fm] 


23 Na 


2.9393 


2.994 


0.523 


127j 


5.5931 


4.749 


0.523 


129 Xc 


5.6315 


4.776 


0.523 


131 Xe 


5.6384 


4.781 


0.523 


132 Xe 


5.6460 


4.787 


0.523 


134 Xc 


5.6539 


4.792 


0.523 


184 w 


6.51(7) 


5.42(7) 


0.535(36) 


186 W 


6.58(3) 


5.40(4) 


0.480(23) 



TABLE VIII: Two-Parameter Fermi coefficients for 23 Na, 127 I, and 129 " 134 Xe as in >3] and for 184 - 186 W as in 20], to be used 
m (ED). <r 2 ) 1/2 , the rms value of the charge radius calculated using the Two-Parameter Fermi distribution with t = 2.30 fm 
for 23 Na, 127 I, and 129_134 Xe, and c and a as given above for 184 ' 186 W, is also given. 



